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We study the possibility of detecting gravitational-waves with matter-wave interferometers, where 
atom beams are split, deflected and recombined totally by standing light waves. Our calculation 
shows that the phase shift is dominated by terms proportional to the time derivative of the grav- 
itational wave amplitude. Taking into account future improvements on current technologies, it is 
^ _ ^ . promising to build a matter-wave interferometer detector with desired sensitivity. 

■ PACS numbers: 04.80.Nn, 04.80.-y, 95.55.Ym 

lO ! I. INTRODUCTION 

(N 

The existence of gravitational vifaves is one of the most important predictions of Einstein's general theory of relativity. 
O The indirect evidence comes from the long term observation of the binary pulsar PSR 1913-1-16 P, Ql- Since 1970s, 
'„ many efforts have been made to find direct evidences. However, it becomes more and more clear that detecting 
5— I ' gravitational waves is really a great challenge to the experimentalists. Two types of detection technologies have been 
developed so far 0, 0]. Resonant mass detectors are the first type of detectors. They were constructed in several 
' laboratories. Some are still running: ALLEGRO in USA, AURIGA and NAUTILUS in Italy, and EXPLORER at 
— ^ CERN. In the meantime, another type of detectors based on laser interferometers are studied and built: LIGO in 
J> . USA, VIRGO in Italy, GEO600 in Germany, and TAMA300 in Japan. These present-day detectors can only reach 
l> ■ the marginal sensitivity to detect gravitational waves from our own galaxy. To reach higher sensitivity and wider 
operating frequency, both types of detectors inevitably have to grow in size, and become more expensive to build. 
Therefore, it is of great worth to study alternative detector designs that can have comparable sensitivity, but with 
smaller size and less expensive. 

■ Actually, the rapid technological advances in atom interferometry provide us a new possibility. Due to their high 
sensitivity, atom interferometers have already been used in many precision measurements, such as the measurement 
of Newton's constant G [H , the test of the equivalence principle Q , and the breakdown of the 1/r^ law at small length 
scales 01 • One can find more details on atom interferometers in many review papers, such as Q. The theoretical 
investigation of detecting gravitational waves with matter- wave interferometry was started by many authors 0-[l2|- 
In recent years, possible experimental schemes based on atom interferometers were extensively studied []"3l - [l9| . These 
schemes use either the super pos ition of two different momentum states of the atom p7l [l8| . or the particle- wave 
duality property of the atom |13l - [l5| . 

In Chiao and Speliotopoulos' scheme [T3| , atom beams are split and recombined by 32nm material gratings, and 
are deflected by crystal mirrors. They called their scheme, MIGO, the matter-wave interferometric gravitational-wave 
observatory. The limitation of their scheme is that only Helium atoms can be used due to the current ability on 
making crystal mirrors. It is preferable to use heavier atoms for detecting gravitational waves. On the other hand, 
it was experimentally known that diffraction gratings of light have more advantages than material gratings. Atom 
interferometers based on standing light waves were successfully constructed a long time ago [13, HH . 

In this paper, we study the possibility of detecting gravitational waves using the matter-wave interferometer based 
on Bragg scattering from standing light waves. Using Bragg scattering, we can split, reflect and recombine a variety 
of atom beams. It is also easier to align standing light waves than material gratings. In Bragg scattering, the atoms 
always remain in a single state so that the atomic phase is not sensitive to external fields and fluctuations in laser 
beams. The rest of the paper is organized as follows. In section II, we calculate the phase shift due to an incoming 
gravitational wave. The difference to laser interferometer detectors can be clearly seen from our result. In section HI, 
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we give a description of our proposed scheme, and analyze main sources of noise that in principle limit the sensitivity 
of our detector. Comments and conclusions will be given in section IV. 

Throughout this paper, we will keep the Planck unit h explicitly, while adopting the usual notation for the speed 
of light c = 1 . Indices are lowered and raised by the Minkowski metric r]^^, . 



II. CALCULATION OF THE PHASE SHIFT 



To make the calculation clear, a schematic diagram of our configuration is depicted in FIG.[T] The length and width 
of the interferometer are denoted hy L// and L± , respectively. Atoms from the beam source with velocity vq are split 
into two beams by angle 2a. The two beams are then reflected, recombined and finally flow into atom detectors. 
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FIG. 1: Schematic diagram of our proposed configuration 

Suppose incoming gravitational waves propagate along the direction orthogonal to the plane of the interferometer, 
which is taken to be the z = plane. In the transverse-traceless gauge [22,] , the metric for this linearized gravitational 
wave can be written as 

ds^ = -dt^ + [S,j + h,j{t - z)]dX'dX\ (1) 

where z = X^, h\ = 0, d^hij = and i^j = 1,2,3. In literature, the metric with hn — — /122 — hsm{2TTft + (po) is 
called the /i+-polarization, and the one with hi2 — /121 — hsm{2'K ft + 0o) is called the /ix -polarization. 

As our measurement is taken in laboratory, it is convenient to work in another coordinate system, related to the 
(t, A')-coordinates by a transformation 

t = t, 

x' = X' + lh]X^. 

2 ■' 

The metric is then rewritten into 

ds^ = -dt^ + Sijdx'dx^ + hij{t - z)[dz - dt\x^dx' + (2) 

where terms of quadratic or higher order in hij have been neglected. In the above transformation, the condition that 
hij(t)x^ and hij are of the same order is assumed. This is necessary for the new coordinates to be consistent with the 
linearized treatment of the gravitational-wave metric. In other words, the condition amounts to require that the full 
size of the gravitational-wave detector should be much smaller than the wavelength of the gravitational wave. 

It is easy to write down the geodesic equations for particles confined in the z = plane in the (t, a;*)-coordinates, 

hy +oihl,v'K,). (3) 



di2 2 

One can check that the physical distances in the z = plane are not affected by gravitational waves. Therefore, this 
coordinate system is called the rigid frame in In this frame, the Lagrangian for a non-relativistic atom becomes 

7T) ... 

L{x\x') = —{x'x^ - h,j£x^ - 2). (4) 
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For a given path, the phase at the final point can be calculated as 



where Si-^j is the action from the initial to the final point, calculated from the Lagrangian Q. 
Before the incoming of gravitational waves, atoms travel along the classical trajectories, 

x,{t) = v,t, yoW = |^|(y_^) forr/2<t<r, 

where — wocosa, T = L///vq, and — iwosina, depending on the path. In the present of gravitational waves, 
both the atom's trajectory and velocity will be perturbed, which can be written as x ~ xq -\- 5x and v = -\- w. To 
solve the geodesic equations for atoms, we need to know the junction conditions at the middle standing light waves. 
The natural conditions are that the geodesic should be continuous and the velocity should satisfy the rigid condition. 
That is, we have the following conditions 

5x\t^T-/2 = Sx\t^T+/2 and w|t_^T-/2 + w\t^T+/2 0. (6) 
Using the same derivation as in [l5| . the phase for each arm is found to be 
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up,down^ = / { J W -\- VqW^ ~hijVQXQ}\up, dowm (7) 

Jo ^0 ^ 



where A^p^ down^ means the phase change in the up and down arms, respectively. Because the /i+-polarization affects 
the up and down arms in the same way, our scheme is only sensitive to the ft. x -polarization. The phase shift is 
obtained by taking the difference of phase changes in the two arms, 

— A<i> = Tr/ZiUgT^ sinatanQC0s(7r/r + (/)o) 

— /iUQrsinatanQ;(l — sin a) sin(7r/r) cos{n fT + 0o) 

+ftworsinatanQ;(l - sina)(l - cos(7r/T)) sin(7r/r + 0o)- (8) 

Note that A$ oscillates in time, because 4>q (the phase of the gravitational wave at the time of the first grating) is 
oscillatory. Let us discuss the result in more detail. Define 

A<f>i — 2nmAf h tana cos{tt fT + (j)o)/{hcosa), 

A<1'2 = —/it'oT' sina tana sin(7r/T + 0o) = 2ft.ij^fcc(B tana sin(7r/r + (/)o), 

where A = i^gT^sinacosa is the area enclosed by the two arms, and kdB = mvo/h is the de Broglie wavenumber. 
In other words, A$i is the contribution involving the derivative of the gravitational wave amplitude, and A$2 is the 
contribution involving the gravitational wave amplitude. We consider three different cases. 

• When /r< 1, 

A$ ^ A$i sin a. 

The phase shift is proportional to A$i, and is very small since a is small. 

• When fT ~ 1, 

A$ A$i + A$2- 
The two different contributions can be comparable. 

• When /r> 1, 

A* A<I>i. 
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The feature of our result is that the phase shift contains terms from two different contributions, which have been 
discussed separately before [Tsj - flGj . In this sense, our derivation is more complete. Especially, in high frequency 
regime 1/7"), the contribution involving the derivative of the gravitational wave amplitude dominates the phase 
shift. 

This phase shift can be measured by detectin g th e oscillation in the intensities of atom beams going into detectors 
A and B. By a similar calculation to the one in [2J], the intensities of atom beams detected by idealized detectors A 
and B are found to be 

r32)-kg+A$)), 

r32)-kg+A$)), (9) 

where r2i = r2 — ri, r32 = — r2, r^'s denote the positions of light gratings, kg is the reciprocal vector associated 
with each light grating, which is assumed to be the same for all light gratings. Its norm is ||kg|| — 27r/a, where a is 
the period of light gratings. Iq is the total intensity of the atom beam from the source. In order to maximize the 
detection sensitivity, we set (r2i — r32) • kg = tt/2. Then it follows that 

Ia - |(l-sin(A$))«|(l-Aci>), 

Ib = |(l + sin(A$))«|(l + Aci>). (10) 

Clearly, the contrast between the atom beam intensities at the two detectors will also help to improve the reliability 
of the detection. 



I A = y(l + cos((r2i - 
Ib = y(l - cos((r2i - 



III. DESIGN CONSIDERATIONS AND NOISE ANALYSIS 



A. Description of our proposed configuration 

The schematic diagram was already shown in FIG.[T] The size of the interferometer is proposed to be L// = 200m 
and Lj_ — Im. The interferometer regions must be in high vacuum about 10~^ Pa in order to keep the coherence of 
atom beams. At this pressure and room temperature, the vacuum contains about n ~ 3 x lO^^m"^ molecules with an 
average velocity of u ~ 500m • s~^. The average collision time between atom beams and air molecules is ~ 80s, 
where the typical value of cross section a is about 10~^®s~^. So the atom beams can keep coherence for about Is. 

The atom beam is emitted from a supersonic source, with velocity vq = 1000m- s~^. It can be well coUimated by 
2D optical molasses. The distribution of the transverse velocity of atoms can be controlled to meet the requirement of 
gravitational-wave detection. Four standing light waves are generated by laser beams from the same laser source. By 
adjusting the intensities of standing light waves, the first and fourth ones can function as 50% — 50% beam splitters, 
and the two in the middle can function as beam mirrors. Interference fringes can be formed in both output atom 
beams of the interferometer. Either beam can be used for detection because they are complementary to each other. 

To yield a good signal-to-noise ratio in our scheme, the supersonic source should have high intensity. For example, 
people have produced Argon beams with high flow rate TZ ^ lO^^atoms/s [24|. To scatter the Ar atom by angle 
2a = 10~^rad, we need a momentum transfer of order lO^hk. Currently, momentum transfer of 24fi,/c has been 
reported (25j . An order-of-10 improvement on large momentum transfer technology is needed. To split, reflect and 
recombine Ar beams of so high flow rate, the estimated laser power in our scheme is of order lO^watts. This high 
laser power can be achieved by using the power amplifier, as in LIGO f26|. 

Taking the above values into cq. ([8]), one can find that A<i> = 2 x 10^ fh. In astrophysics, there exist many 
interesting sources of gravitational waves, such as compact binary coalescences (CBCs) of neutron stars or black holes 
[13]. The gravitational- wave emitting from CBCs will last for over tens of seconds with frequency in the 10 — lO'^Hz 
range, and amplitude h < 10^^^. Our detector would produce a phase shift A"I> ^ 10^^^ — 10~^rad. Whether this 
phase signal can be detected depends on various sources of noise in our scheme, which can have more or less effects 
on the detection ability of the detector. A thorough analysis of noise is necessary for the detector design, which will 
be reported in our future work. Here, we focus on shot noise and seismic noise. 
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FIG. 2: Amplitude spectral density curve for shot noise 



B. Shot Noise 



Shot noise is the fundamental limitation on the sensitivity of any atom interferometer. The fluctuation in the phase 
shift measurement is inversely proportional to the square root of the total number of atoms, 
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(11) 



where TVq is the total number of atoms, and ta is the experimental time. This determines the minimal detectable 
gravitational-vifave signal of the detector. To make contact with other literature on gravitational-wave detection, it 
is useful to write down the corresponding amplitude spectral density, which is independent of the experimental time. 
The formula is found to be 



where 



''shot 



if) 



2h 



m\Cif)\VTl 



(12) 



C{f) = [-7r-(;o7'^sinatanasin(7r/T)/ + forsinatana(l - sina) sin^(7r/r) 
+iiQTsinatana(l - sina)(l - cos(7r/T)) cos(7r/r)] 

+i[TrvQT'^ sinatanacos(7r/T)/ — UQTsinatana(l — sina) sin(7r/T) cos{tt fT) 
-|-i'QTsinatana(l — sina)(l — cos(7r/r)) sin(7r/T)]. 



(13) 



Using the proposed scheme parameters, the amplitude spectral density curve is drawn in FIGl2] The curve clearly 
shows that the sensitivity of our scheme gets better at higher frequencies. 
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C. Seismic Noise 

Seismic noise is another important factor that will seriously affect the detection ability of the gravitational-wave 
detector. Typical seismic noise levels are psj 



^ seisvmc 



if) 



lO-^m/VHz 
{lOHz/ff ■ 10"Sii/\/Hz 



for \Hz < f < lOHz 
for / > lOHz 



(14) 



In our scheme, the seismic noise has two main effects on the detection. One effect is that it will change the 
phase shift A$. However, this effect is very tiny, for example, at the frequency of lOOHz, the seismic noise is about 
10^^^m/\/Hz. The interferometer area will change about AA = 10~^m^/\/Hz by this seismic noise. Then the phase 
shift will change about 10~^°rad for ta = 10s, which is really small. 

The other effect is that the seismic noise will break the condition, (r2i — r32) • — 7t/2, we assumed before. If we 
don't adopt any vibration isolation system, the corresponding amplitude spectral density is 
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(15) 



Again, using the above scheme parameters and a = A/2 ^ 400nm, we can draw the amplitude spectral density 
curve, as shown in FIG (3] Seismic noise should be reduced to be smaller than shot noise, which means 



A(r2i - r32) • kg < 



(16) 



Compare FIG|3]with FIG 121 we find that we need an isolation factor of about 10^ in the 10 — lOOOHz range, which 
has already been achieved in LIGO [26| . 
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FIG. 3: Amplitude spectral density curve for seismic noise 
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IV. CONCLUSIONS 



Through this prehniinary study, we find that in high frequency regime the phase shift in our proposed scheme is 
proportional to the time derivative of the gravitational wave amplitude. This is a good feature of our scheme, which 
makes it suitable for detecting gravitational waves in the frequency band 10 — lOOOHz, in which there exist a lot of 
interesting astrophysical sources. Considered current technologies and their future improvements, it is possible for 
our proposed detector to reach a high sensitivity at a much lower cost. For all these good features, our scheme is a 
good candidate scheme for building future gravitational-wave detectors. 
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